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ABSTRACT 

Despite the great potential of the [/-band galaxy luminosity function (GLF) 
to constrain the history of star formation in clusters, to clarify the question of 
variations of the GLF across filter bands, to provide a baseline for comparisons 
to high-redshift studies of the cluster GLF, and to estimate the contribution of 
bound systems of galaxies to the extragalactic near-UV background, determi- 
nations have so far been hampered by the generally low efficiency of detectors 
in the [/-band and by the difficulty of constructing both deep and wide sur- 
veys. In this paper, we present [/-band GLFs of three nearby, rich clusters to 
a limit of M\j « —17.5 (Mjj + 2). Our analysis is based on a combination of 
separate spectroscopic and -R-band and [/-band photometric surveys. For this 
purpose, we have developed a new maximum-likelihood algorithm for calculating 
the luminosity function that is particularly useful for reconstructing the galaxy 
distribution function in multi-dimensional spaces (e.g., the number of galaxies as 
a simultaneous function of luminosity in different filter bands, surface brightness, 
star formation rate, morphology, etc.), because it requires no prior assumptions 
as to the shape of the distribution function. 

The composite luminosity function can be described by a Schechter func- 
tion with characteristic magnitude Mjj = —19.82 ± 0.27 and faint end slope 
au = —1.09 ± 0.18. The total [/-band GLF is slightly steeper than the i?-band 
GLF, indicating that cluster galaxies are bluer at fainter magnitudes. Quiescent 
galaxies dominate the cumulative [/-band flux for Mjj < —14. The contribution 
of galaxies in nearby clusters to the [/-band extragalactic background is < 1% 
Gyr -1 for clusters of masses ~ 3 x 10 14 to 2 x 10 15 M Q . 

Subject headings: galaxies: clusters: general — galaxies: evolution — galaxies: 
luminosity function, mass function — methods: statistical 



1 Steward Observatory, The University of Arizona, 933 N Cherry Ave, Tucson 85721 AZ 
2 Astronomy Department, University of Massachusetts, Amherst, MA 01003 



- 2 - 



1. INTRODUCTION 

The galaxy luminosity function (GLF) is one of the most fundamental statistics of 
galaxy populations. Its shape and variation with environment provide a crucial constraint 
on any model of galaxy evolution. Recent studies of the GLF in clusters of galaxies, based 
on spectroscopic surveys, have been carried out in the i?-band (Christlein & Zabludoff 2003) 
and the frj-band (de Propris et al. 2003). However, there is a large amount of additional 
information to be gleaned from [/-band GLFs of cluster galaxies. 

Why do existing luminosity functions in the R- and ftj-band not provide us with a 
complete picture of the galaxy population in clusters? Different filter bands are sensitive to 
different stellar populations. For example, blue and near-UV bands are most sensitive to 
recent star formation, while red and near-IR bands better approximate the total stellar mass. 
Determinations of the GLF in different magnitude bands are therefore complementary, and 
the [/-band in particular promises insight into a number of important questions: 

First, several studies (Bromley et al. 1998; Madgwick et al. 2002; Christlein & Zabludoff 
2003) show that star-forming and quiescent galaxies have very different GLFs. It is therefore 
reasonable to expect that star formation would affect the shape of the overall GLF, and that 
it would do so differently in different filter bands. Other phenomena, such as the presence of 
dust or the metallicity of a galaxy population, also affect galaxy colors and could introduce 
inhomogeneities in comparisons between different filter bands. Understanding how strongly 
such color-dependent effects influence the GLF is crucial to evaluating the significance of 
discrepancies between observational determinations of the GLF in different filter bands, and 
between observed GLFs and model predictions. 

Second, with the availability of large telescopes and new detection techniques (Madau 
et al. 1996), recent years have seen growing interest in observations of clusters at high 
redshift. Surveys of the GLF in such systems reach to approximately 2 mag fainter than the 
characteristic bright magnitude M* at redshifts up to z « 1 (Stanford et al. 2002; Massarotti 
et al. 2003). High- redshift observations in red and near-IR bands typically observe blue or 
near-UV rest-frame wavelengths (Fig. 1). Measuring the evolution of such fundamental 
statistics as the GLF in these bands has been complicated by the lack of corresponding low- 
redshift GLFs. Our [/-band GLF for clusters at z ~ provides a baseline for comparisons 
to high-redshift clusters that are observed in the [/-band rest-frame band with a comparable 
depth (< M* + 2). 

Third, there is an extragalactic background in the [/-band at a level of 2 to 4xl0 -9 erg 
s _1 cm~ 2 sr _1 A -1 (Bernstein, Freedman & Madore 2002; Henry 1999). It is thought that a 
majority of this light is produced by stellar nucleosynthesis rather than AGN or other non- 
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Fig. 1. — Comparison between rest-frame UBVRI and observed UoBoVoRolo filter bands at 
different redshifts. Tilted lines show how the rest-frame bands are mapped onto observed 
bands (vertical solid lines) for any given redshift. At z=0.8, observations in the i?-band 
detect the rest-frame [/-band. 
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stellar sources. Therefore, the [/-band extragalactic background is of cosmological interest 
as a cumulative constraint on stellar nucleosynthesis. Furthermore, it is related to the far- 
UV background, knowledge of which is crucial to understanding the ionization state of the 
intergalactic medium. From galaxy number counts in deep images, numerical estimates of the 
contribution of resolved and unresolved normal galaxies to this background have been made 
(Bernstein, Freedman & Madore 2002b). Our present work makes it possible to calibrate a 
relation between cluster mass and cumulative [/-band luminosity. By coupling this with a 
cluster mass function (Jenkins et al. 2001), we can make an independent estimate for the 
contribution of nearby clusters (with a mass range of 10 14 to 10 15 M Q ) to the extragalactic U- 
band background light. Furthermore, assuming that there is no break in the mass-luminosity 
relation, we can extrapolate it to roughly estimate the contribution from lower-mass, group- 
like systems (10 12 to 10 14 M ). 

Determinations of the [/-band GLF in clusters at low redshift have so far been com- 
plicated by the generally low efficiency of detectors in the [/-band and by the challenges of 
surveying both wide and deep. Beijersbergen et al. (2002) give a deep determination in the 
Coma cluster, but use statistical background subtraction rather than a spectroscopically se- 
lected sample to account for contamination of the sample by fore- and background galaxies. 
As Valotto, Moore & Lambas (2001) have demonstrated, this technique can be subject to 
large systematic errors if the background is inhomogeneous. Spectroscopic samples allow for 
superior background subtraction. 

The recent availability of blue-sensitive wide-field detectors now makes measurement 
of [/-band GLFs possible. The combination of a [/-band photometric survey (Mcintosh, in 
prep.) with a spectroscopic sample of cluster galaxies from the Las Campanas Nearby Cluster 
Survey (Christlein & Zabludoff 2003) enables us to, for the first time, present [/-band GLFs 
of three nearby clusters from a spectroscopically selected galaxy sample. The availability 
of i?-band photometric data for the same set of galaxies allows us to make self-consistent 
comparisons between i?-band and [/-band GLFs from the same sample. 

Our procedure of combining additional photometric data in the [/-band with an existing 
survey whose completeness is known in the i?-band requires us to deal with at least three 
quantities in the calculation of the GLF: the apparent [/-band magnitude {mjj), the i?-band 
magnitude (m^), and the i?-band surface brightness (/i_r). We have therefore developed 
a new algorithm for the calculation of the GLF that is particularly suited for such multi- 
variate analyses. This new algorithm is a variant of maximum likelihood estimators and 
retains their advantage of being unbiased by density inhomogeneities due to large scale 
structure. In addition, it offers the benefit that no analytical or binned form for the galaxy 
parent distribution needs to be assumed a priori, a great advantage in multi-variate problems, 
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where analytical forms often do not exist and binning the distribution in a multi-dimensional 
space is inefficient. 

We discuss our cluster sample and the spectroscopic and photometric surveys in §2. In 
§3 and App. A, we introduce our new GLF algorithm, the Discrete Maximum Likelihood 
method, and we discuss the completeness of each survey, as well as systematic corrections 
to account for biases related to color terms in the sampling fraction in App. B. §4 gives 
our results and discussion. We present the GLFs for the three individual clusters in §4.1. 
We determine the composite GLF for all galaxies as well as for emission line (star forming 
and active) and non-emission line (quiescent) galaxies in §4.2. We then compare the U- and 
i?-band GLFs, which are calculated from the same sample and from the same processing 
pipeline (§4.3). In §4.4, we examine the contribution of clusters to the metagalactic [/-band 
background. In Appendices C and D, we discuss the effects of the spatial and magnitude 
limits of our survey on our ability to sample the total [/-band flux from clusters. Our 
conclusions are presented in §5. 

2. THE DATA 

2.1. The Cluster Sample 

Our sample consists of three clusters (Abell 496, Abell 754, Abell 85) from the original 
spectroscopic and i?-band imaging sample of Christlein & Zabludoff (2003). The Christlein 
& Zabludoff clusters were selected based on 1) their visibility from Las Campanas, 2) the 
availability of some prior spectroscopic and X-ray data in the literature, 3) their redshifts, 
which allowed us to sample the cluster to at least one virial radius with the 1.5° xl.5° field 
of the fiber spectrograph field, and 4) their range of velocity dispersions, which suggest a 
wide range of virial masses. The properties of these clusters (for H = 71 km s _1 Mpc -1 , 
VL m = 0.27 and il\ = 0.73, as applied throughout this paper) are given in Table 1. We refer 
readers to Christlein & Zabludoff (2003) for details on the spectroscopic survey and data 
reduction. 



2.2. i?-band Survey 

We derive our master galaxy catalog from a photometric survey of the clusters in the 
i?-band (Christlein & Zabludoff 2003). This catalog is complete within certain magnitude 
and surface brightness limits, so we use it as the reference to estimate the completeness of 
the spectroscopic and the [/-band photometric catalogs as a function of (mR,iiR). Details 
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of this survey, the image reduction, photometry, and the construction of the master catalog 
are in Christlein & Zabludoff (2003). 



The f/-band cluster galaxy data come from wide-field (1° xl° ) imaging of A85, A496, 
and A754 during a January 2000 run using the NOAO Mosaic Imager on the Kitt Peak Na- 
tional Observatory (KPNO) 0.9-meter Telescope. Complete details of the sample selection, 
observations, reductions and photometric calibrations are reported in Mcintosh, Rix, Cald- 
well & Zabludoff (2004, in prep.). The final sample contains [/-band and V^-band magnitudes 
for a total of 631 spectroscopically confirmed cluster member galaxies, and is comparable in 
depth and in membership coverage to the most comprehensive study of the Coma cluster 
by Terlevich, Caldwell, & Bower (2001). In physical units, the [/-band survey region covers 
about 1.7 x 1.6 Mpc 2 in A496, 2.7 x 2.6 Mpc 2 in A754, and 2.7 x 2.6 Mpc 2 in A85. 

The data are well-flattened and carefully corrected to ensure uniformity of the photomet- 
ric zero point using a customized reduction pipeline that follows standard image reduction 
techniques and uses the IRAF 1 environment. We perform the basic reduction of the [/-band 
data using the IRAF mscred package. The images are spatially and spectrally flattened 
(< 1% deviations globally) using an optimized night-sky flat-field frame constructed from 
individual science exposures with all objects masked. We calibrate the data astrometrically 
to the USNO-v2.0 system (Monet et al. 1996) and photometrically to Johnson (1966) U 
magnitudes on the Landolt (1992) system. 

We perform source detection and extraction using SExtractor (Bertin & Arnouts 1996) 



1 IRAF is distributed by the National Optical Astronomical Observatories, which are operated by AURA, 
Inc. under contract to the NSF. 



2.3. [/-band survey 



Table 1. The Cluster Sample 



Cluster 



N 



cz [km/s] 



Am [mag] cz range [km/s] a [km/s] 



center (J2000.) 



area arcmin 



A496 
A754 
A85 



241 
415 
280 



9910 ± 48 
16369 ± 47 
16607 ± 60 



35.78 
36.90 
36.94 



7731 - 11728 728 ±36 04:33:37.84,-13:15:44.5 
13362 - 18942 953 ±40 09 : 08 : 32.00, -09 : 37 : 00.0 
13423 - 19737 993 ±53 00 : 41 : 50.46, -09 : 18 : 11.6 



59.9 x 58.7 

60.1 x 58.5 

60.2 x 58.3 



Notes: N is the number of sampled galaxies per cluster, cz is the mean rcdshift, Am the distance modulus. Center 
coordinates and survey area are for the [/-band photometric survey, which has a smaller coverage than the i?-band 

photometric or the spectroscopic survey. 
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with the following configuration parameters defining our definition of an imaged source: 
minimum of 5 detected pixels (DEBLEND_MINAREA) above a background threshold of 
3<Tbk g (DETECT_THRESH), with overlapping sources deblended if the contrast between flux 
peaks associated with each object is > 0.05 (DEBLEND JV11NCONT) . We confirm that these 
parameters provide good source detection and deblending by visually inspecting random 
regions from each image. We reject sources flagged (FLAGS> 4) as saturated or otherwise 
bad, and we exclude all sources within 1 arcmin of image edges. The empirical magnitude 
limits where the £7-band source count distributions flatten and turn over are roughly 20.7 
(A85), 21.0 (A496), and 20.5 (A754) mag. We separate stellar and extended sources by 
fitting a PSF-convolved bulge+disk model to the light profiles of detected sources using 
G1M2D (Simard et al. 2002) following the method described in Mcintosh et al. (in prep.). 
This method is robust to V = 19 mag, which corresponds to ~ 20.5 in U for most cluster 
members. 

From 975 redshift coordinates within the coverage of the three cluster U images, we find 
727 cluster members (with czi = (cz) cl ± 3er c i) and 248 with recessional velocities outside 
of these bounds. We cross-correlate the coordinates of member galaxies from the redshift 
data with the [/-band source positions and achieve 631 members from our imaging catalog. 
We define image/redshift matches to be the nearest within 5" and find the mean coordinate 
separation is < 2". Only 2.4% (12) of the 88 redshift positions without a U image detection 
are brighter than U = 19.5. Finally, 112 U < 19.5 extended sources have no redshifts due to 
the fractional incompleteness of the spectroscopic sample at R > 16 (Christlein & Zabludoff 
2003). Therefore, our [/-band galaxy catalog from clusters A85, A496, and A754 is ~ 85% 
complete down to U ~ 19.5 mag. 

We use the dust maps by Schlegel, Finkbeiner & Davis (1998) for extinction corrections 
and apply a k-correction of K v = 0.065(^/0.02), as appropriate for early-type galaxies 
(Pence 1976), to the [/-band magnitudes. We apply this correction to all galaxies in this 
sample, because cluster galaxies are predominantly early types, individual morphological 
classifications are not available, and the k-correction for extreme late types differs from the 
value for early types by less than 0.1 mag even for the most distant sample cluster. 

3. CALCULATING THE LUMINOSITY FUNCTION 

3.1. A New Method 

For reconstructing the luminosity function of our survey, we use a new statistical method, 
which we refer to as the Discrete Maximum Likelihood Method. The DML retains the 
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advantages of Maximum Likelihood estimators in being unbiased by density inhomogeneities 
in a sample and is therefore easily applicable to field as well as cluster samples. In contrast to 
most maximum likelihood-based LF algorithms (Sandage, Tamman & Yahil 1979; Efstathiou, 
Ellis & Peterson 1988; Blanton et al. 2001), which which were developed when luminosity 
alone was the primary variable of interest in determining galaxy distribution function, the 
DML algorithm does not assume an ansatz (binned or analytic) for the distribution function 
a priori, and is therefore independent of the dimensionality of the parameter space of interest. 
It is therefore ideally applicable to multivariate distribution functions, i.e., the abundance of 
galaxies as a simultaneous function of luminosity in different filter bands, surface brightness, 
environment, star formation, morphology, etc. The derivation of the Discrete Maximum 
Likelihood Method is discussed in detail in Appendix A. 



3.2. The Sampling Fraction 

The reconstruction of the GLF, as described in App. A, requires knowledge of the 
sample completeness, f(xh \ Fi). f is the probability that we have both [/-band photometry 
and spectroscopic information for a given object with certain physical properties x n (e.g., 
absolute magnitude and surface brightness) if it is in a particular field Fi (characterized 
by redshift, Galactic foreground extinction, position on the sky). Due to the design of 
our survey, the sampling fraction is not known analytically, but has to be reconstructed 
from the data. Our i?-band detection catalog serves as the reference against which we 
calculate the completeness of the spectroscopic and [/-band photometric samples. Appendix 
B describes the determination of the sampling fraction, as well as systematic corrections that 
are necessary to account for color selection effects. 



3.3. Consistency checks 

The introduction of a new method for the calculation of the luminosity function requires 
us to demonstrate its consistency with established algorithms. We confirm this by comparing 
the composite i?-band GLF obtained from the SWML method (as used in Christlein & 
Zabludoff (2003)) with our new algorithm. Fig. 2 shows these two GLFs (upper curves, 
unnormalized) , determined for the three clusters in our survey. The GLFs are statistically 
indistinguishable. The results of both methods are equivalent for calculating distribution 
functions in low-dimensional parameter spaces, but the DML offers greater efficiency and 
convenience in treating multi-dimensional problems. 
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Fig. 2. — Upper pair of curves: Comparison between composite i?-band cluster GLFs using 
the Stepwise Maximum Likelihood method (dashed line; Efstathiou, Ellis & Peterson 1988) 
and our Discrete Maximum Likelihood method (solid line), determined over the area of the 
i?-band photometric sample. Normalization is arbitrary. Lower pair of curves: Comparison 
between GLFs from all galaxies within [/-band survey area (dashed line) and from those 
galaxies with [/-band photometry only, after completeness correction (solid line). Normal- 
ization is galaxies / magnitude / cluster. 
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We also test whether our sampling fraction model accounts correctly for the incomplete- 
ness of the [/-band photometric sample. This is particularly important because we model 
the completeness of the {/-band catalog as a function of (m R ,ii R ), while the strongest de- 
pendence should be on (mu,fj,u). We calculate the composite i?-band GLF from all galaxies 
in the i?-band master catalog that lie within the spatial boundaries of the {/-band survey, 
regardless of whether we know their {/-band magnitude. We then calculate an i?-band GLF 
from galaxies for which we have U -band photometry, and apply the appropriate completeness 
corrections. The agreement between the two GLFs is again excellent, indicating that the 
incompleteness of the {/-band catalog does not introduce systematic biases. We emphasize 
that this i?-band GLF is plotted using the same weighting and normalization factors as the 
{/-band GLF that we present below. 



3.4. Coverage of (M R ,M V , Hr) Parameter Space 

Our survey has magnitude limits in two different filter bands, R and U. At each Mu, 
the R band magnitude limit introduces an incompleteness in the calculated {/-band GLF, 
because, for any given Mu, some fraction of galaxies are too faint in M R to have been 
sampled. To avoid this problem, we calculate {/-band GLFs only down to a limit at which 
the fraction of lost galaxies is likely to be small. 

We use Fig. 3 to identify this completeness limit for our {/-band GLF. The figure shows 
the completeness-corrected galaxy distribution (greyscale) in the (Mu, M R ) plane. Dots 
represent the individual galaxies in the sample. The bold solid line shows the faintest M R 
to which a galaxy with a given Mu and any /j, R typical for galaxies at that Mu could have 
been sampled. The effective absolute magnitude limit at the faint end is M R « —17.3. Fig. 
3 shows that, for Mu > —18, the galaxy distribution has not been sampled completely; in 
any Mu bin substantially fainter than that limit, a substantial fraction of galaxies are likely 
to be missing due to the M R limit. We adopt Mu = —17.5 as the absolute magnitude limit 
for our analysis of the {/-band GLF; at this limit, most of the galaxy distribution function 
has been sampled, and only a small fraction of galaxies (~ 5%, if the galaxy distribution 
in M R for a given Mu is approximately Gaussian) in the faint-M^ tail are likely to be lost. 
The two more distant clusters, A754 and A85, have only been sampled to M R = —19.0. 
Therefore, we consider the individual {/-band GLFs for these clusters to be reliable down to 
Mu = -18.5. 

The projection of the distribution onto the (Mu, M R ) plane is noticeably broader at 
fainter magnitudes than at the bright end, a result of the greater fraction of blue, emission 
line galaxies. 
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Fig. 3. — Galaxy distribution in the (Mjj, Mr) plane in greyscale as calculated using the 
Discrete Maximum Likelihood method. Points denote cluster members with both redshifts 
and [/-band photometry available. The bold horizontal line at Mr m —18.0 gives the effective 
absolute magnitude limit in R, i.e., the Mr at which the sampling fraction is zero for a galaxy 
of a given My and typical fiR at that M\j. Galaxies fainter than M\j m —17.5 are too faint 
in Mr to have been sampled, thus motivating our absolute magnitude limit of M v = —17.5. 
The more distant clusters are only sampled to Mr —19.0 (thin horizontal line), motivating 
a limit of M v = -18.5. 
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4. RESULTS AND DISCUSSION 

4.1. [/-band GLFs for Individual Clusters 

Fig. 4 shows the individual GLFs for the three clusters (calculated by the method in 
§3.1) superimposed on the composite GLF, which we calculate in §4.2. Table 2 gives the 
numerical values for the cluster GLFs down to the effective sample limits. 

The shapes of the individual cluster GLFs are consistent with each other within la 
under a x 2 test if the normalization is adjusted to optimize the agreement. However, the 
cluster GLFs are significantly offset in normalization. By drawing 1000 Monte Carlo samples 
of galaxies from the three cluster GLFs, we find the expected number of galaxies at M v < 
-18.5: 91.8±11.4 in A496, 222.0±17.4 in A754, and 143.1±13.8 in A85 over their respective 
survey regions. If we truncate the survey regions to comparable physical radii (taking into 
account both the angular diameter distance and the fact that characteristic length scales 
tend to increase linearly with a (Girardi et al. 1998)), the numbers are 91.8 ± 11.4 in A496, 
189.2 ± 15.8 in A754, and 120.9 ± 12.9 in A85. The differences between A754 and the other 
two clusters are significant and most likely reflect the higher mass of A754. 

In the next section, we will form a composite GLF from the three individual clusters, de 
Propris et al. (2003) find evidence that the shape of the GLF in clusters varies as a function 
of distance from the cluster center, raising the question whether it is legitimate to form 
a composite GLF from three clusters with inhomogeneous sampling radii. In our sample, 
we find no statistically significant differences in the shapes of the GLFs of A754 and A85 
between the truncated and untruncated samples, which legitimizes our approach to forming 
the composite GLF. 

4.2. Composite [/-band GLFs 

We apply the DML algorithm to the complete sample of three clusters and then bin the 
recovered galaxy distribution function over M v . We repeat this process for subsamples of 
emission line (EL) and non-emission line (NEL) galaxies separately. EL galaxies are defined 
as having an equivalent width of the [Oil] A3727 doublet of 5A or more and are thus star- 
forming or active galaxies. NEL galaxies have an equivalent width of less than 5A and are 
hereafter classed as quiescent galaxies. Fig. 5 shows all three composite £/-band GLFs. We 
have applied corrections for B-R color terms, as described in App. B, in all cases, la error 
bars are based on Poisson errors. Table 3 gives the numerical values for the three GLFs, 
including the numbers of galaxies in each bin. All GLFs are normalized as described in App. 
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Fig. 4. — [/-band GLFs for individual clusters. Normalization of all GLFs is preserved. We 
show the GLFs for A496 (open squares), A754 (crosses), and A85 (triangles). Solid line 
shows best fit Schechter function to composite GLF. 
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A to represent the average number of galaxies per magnitude and cluster within the region 
of the ?7-band photometric survey. 

We have fitted Schechter functions (Schechter 1976) to these GLFs. The functional form 

is 

0(M) = 0*(lO a4 ( M *- M )) 1+a ea;p(-lO a4 ( M *- M )), (1) 

where the nomenclature is used to distinguish this function, which is only a function of 
absolute magnitude, from the (f in Eqn. Al, which is a function of an argument of arbitrary 
dimensionality. 

The faint end slope is not well-constrained in this sample because the GLFs are only 
sampled to a fairly bright magnitude limit. The best fit Schechter parameters are M v = 
-19.821°;^, au = -1.09+°;}^, and = 142 ±5 for the overall GLF. Errors in $y are for the 
best fit a and M* held fixed. Corrections for the B-R color terms in the sampling fraction, 
as described in App. B, are << 1% in and thus negligible for the overall GLF. 

The NEL GLF is also described well by a Schechter function. The best fit parameters 
are M* = -19.77l{j$j, a v = -0.97t°;^ and <f>* v = 133 ±6 for the NEL GLF. Corrections for 
the B-R color terms are more important here, because the color difference between EL and 
NEL galaxies is larger than between field and cluster galaxies. Without our color corrections, 
the faint end slope would be shallower by Aa 0.11. 

For the EL galaxies, both a Schechter function and a power law provide acceptable fits. 
To reduce the number of free parameters, and because the Schechter parameters are very 
weakly constrained by the EL GLF, we decide to fit a power law of the functional form 

0(M) = 0*(lO°- 4 (-(M+2O)) ) l+ a (2) 

The best fit parameters are au = — 2.16lo!i9 and 0* = 4.45+q;52- Again, color corrections are 
important here: the uncorrected EL GLF is steeper by Aa m 0.12. 

4.3. Comparison to i?-band GLFs 

Earlier attempts to determine the GLF in clusters have yielded conflicting results, with 
reported faint end slopes varying from a ~ — 1.0 to a ~ —2.2. One possible explanation of 
these discrepancies could lie in the use of different filter bands. Comparisons between two 
recent, spectroscopically selected studies of the cluster GLF (Christlein & Zabludoff 2003; 
de Propris et al. 2003) do not show evidence for systematic variations in the faint end slope 
a between the i?-band and the b j-band. Our present work provides an even more stringent 
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Table 2. [/-band GLFs. Corrected number of galaxies per mag and cluster. Number of 

sampled galaxies per bin in parentheses. 



Mu 


A496 


A754 


A85 


-22.25 




2.04(1) 


1.98(1) 


-21.75 


2.02(1) 






-21.25 




6.21(3) 


3.98(2) 


-20.75 


4.03(2) 


19.28(9) 


12.38(6) 


-20.25 


18.31(8) 


42.88(21) 


38.42(18) 


-19.75 


37.28(17) 


70.49(34) 


55.15(26) 


-19.25 


66.42(30) 


122.43(56) 


90.64(42) 


-18.75 


55.13(22) 


180.50(73) 


82.94(36) 


-18.25 


57.23(24) 


145.09(50) 


92.78(38) 


-17.75 


73.61(25) 


177.61(44) 


111.43(36) 



Table 3. Composite [/-band GLFs. Corrected number of galaxies per mag and cluster. 
Number of sampled galaxies per bin in parentheses. 



Mu 


all 


EL 


NEL 


-22.25 


1.29(2) 


0.62(1) 


0.66(1) 


-21.75 


0.64(1) 


0.63(1) 




-21.25 


3.22(5) 




3.31(5) 


-20.75 


11.22(17) 


1.33(2) 


10.29(15) 


-20.25 


31.53(47) 


2.40(3) 


29.78(43) 


-19.75 


52.03(77) 


4.98(7) 


48.11(68) 


-19.25 


90.73(128) 


15.65(19) 


75.71(104) 


-18.75 


112.20(131) 


26.65(22) 


83.11(100) 


-18.25 


110.55(112) 


26.15(12) 


96.40(89) 


-17.75 


161.86(105) 


37.70(9) 


110.45(70) 
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Fig. 5. — [/-band GLFs for all galaxies (circles), non-emission line galaxies (crosses) and 
emission line galaxies (open squares). Data points for emission and non-emission line galaxies 
are displaced by ±0.05 mag for display purposes only. Schechter functions for overall (bold) 
and non-emission line galaxies, as well as the power law fit for emission line galaxies, are 
also shown. 
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test for the uniformity of the GLF across different filter bands, because the [/-band is even 
more sensitive to star formation and dust than the bj band. If differences in the GLFs exist, 
they are likely to be revealed in a comparison between U and R. 

We first compare the [/-band and i?-band GLFs using a \ 2 test, allowing for one degree 
of freedom (the shift in M) to minimize x 2 . This approach not only tests the GLFs for 
consistency in shape, but also provides a numerical estimate of the magnitude offsets between 
the ?7-band and the i?-band. The magnitude offset that optimizes \ 2 is U — R = 1.56 ±0.04. 
The maximum probability of consistency is p = 0.10. Although the R- and [/-band GLFs 
are formally consistent under this \ 2 test, a comparison of the ratio of galaxies with M v < 
My < —17.5 to galaxies with My < M v reveals that the [/-band GLF is slightly steeper; the 
ratio is 8.3 ±0.9 for the [/-band versus 5.2 ±0.5 for the i?-band. Fig. 6 shows a superposition 
of the [/-band GLF and the shifted i?-band GLF. 

The difference in the shape of the two GLFs indicates that there is a color gradient, 
d %~R^ ' * n ^ ne sample: fainter cluster galaxies are bluer. What is the cause of this color 
gradient? Factors that could influence (U — R) colors are star formation, nuclear activity, 
dust, and metallicity. Of these, AGN do not appear to play a significant role in our study: 
the NASA Extragalactic Database lists only one Seyfert galaxy in our sample. We now test 
whether differences in star formation activity alone could account for the color gradient. For 
this purpose, we separate our sample into star-forming (EL) and quiescent (NEL) galaxies. 
If star formation activity is solely responsible for the color gradient, then allowing for star- 
forming and quiescent galaxies to have different colors should account for the discrepancy in 
the overall GLF. 

Fig. 7 compares the U- band and i?-band GLFs of the EL and NEL subsamples. The 
optimal magnitude shift that minimizes \ 2 f° r t ne EL GLFs is 0.76 ± 0.30. The maximum 
X 2 probability is 0.89, indicating that these GLFs are consistent with each other in shape. 
For the NEL galaxies, the \ 2 comparison yields a magnitude shift of U — R = 1.65 ± 0.13 
with a probability of p — 0.97. The ratio of faint to bright galaxies, determined as above 
and relative to the same absolute magnitude, M V) is not significantly different between the 
U- and the i?-band in either the EL or NEL subsample, suggesting that U — R colors are 
more homogeneous in the subsamples than in the overall sample. 

If we shift each of the i?-band EL and NEL GLFs by their respective best U — R magni- 
tude offset and combine them, the resulting GLF has a steeper faint-to-bright ratio than the 
original i?-band GLF, indicating that star formation accounts for some of the discrepancy 
between the U- and i?-band GLFs. However, the shifted i?-band and the [/-band GLF are 
still significantly different with regard to the ratio of faint to bright galaxies. Therefore, our 
simple bisection of the sample into star-forming and quiescent galaxies by [Oil] EW does not 
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Fig. 6. — Superposition of [/-band and i?-band GLFs. The i?-band GLF (dotted line) has 
been shifted by A = U — R = 1.56 mag to maximize the agreement between both GLFs. 
Detailed analysis shows the ?7-band GLF to be steeper than the i?-band GLF. 
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Fig. 7. — Same as Fig. 6 for emission line (top) and non-emission line (bottom) galaxies. 
The i?-band emission line GLF has been shifted by A = U — R = 0.76, and the non-emission 
line GLF by A = U — R = 1.65. Within each subsample, we find no significant discrepancy 
in the GLF shapes between the U- and i?-band. 
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completely explain the color gradient. Other effects, such as dust and metallicity gradients 
with luminosity, are probably responsible for the residual color gradient, or else the above 
binning is not fine enough to homogenize colors in each subsample. Clarification of this 
question will have to wait for future investigations with larger surveys. 

Color gradients can introduce differences in the faint end slope between different filter 
bands. Can these differences account for the wide range of a that has been reported in the 
literature? To answer this question, we simultaneously fit Schechter functions to both the 
U- and i?-band GLFs to determine the difference ajj — cxr and its uncertainty interval. If 
we fix My — M* R at its best fit value of 1.93, we obtain ajj — (Xr = +0.03lon- This is 
consistent with the expectation that color gradients within the galaxy population introduce 
only a small difference of a v — a R = — (ocr + 1)- With an estimated ^^"^ = —0.08, 

based on the U — R color offsets for EL and NEL galaxies above and the EL/NEL ratios from 
Christlein & Zabludoff (2003), the expectation value for our sample is au — ctR = —0.01. 
Therefore, the faint end slope a in clusters is only a weak function of the filter band. This 
is qualitatively consistent with Paolillo et al. (2001), who found the faint end slope a to be 
nearly identical in the i-, r-, and g-bands. 



4.4. Contribution from Clusters and Groups to Near-UV Background 

Understanding the contribution from normal galaxies to this background light provides 
important constraints on the star formation history of the universe and the relative impor- 
tance of non-stellar sources, such as AGN or Ly-a recombination radiation from intergalactic 
gas (Tyson 1995). Observational estimates for the extragalactic background intensity in the 
near-UV are on the order of 2 — 4 x 10~ 9 erg s^ 1 cm~ 2 sr _1 A -1 (Bernstein, Freedman 
& Madore 2002; Henry 1999). Past studies (Pozzetti et al. 1998; Bernstein, Freedman & 
Madore 2002b) have estimated the contributions from field galaxies to this background ra- 
diation from number counts in the Hubble Deep Field. Our data provide an independent 
estimate of the contribution of nearby massive clusters, which are different from the field 
with regard to their constituent galaxy populations and evolutionary history. 

What is the contribution of clusters such as those in our sample to the extragalactic 
background in the [/-band? The rate of increase of the [/-band background intensity from a 
population of sources with differential luminosity L UjX and differential spatial density dn/dM 
(where M is the mass) is 

dl\ c f r dn , 

IT = S J L "^lM dM - < 3) 

An estimate of the mass function, dn/dM (M), is given by Jenkins et al. (2001) from 
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numerical simulations. We use their mass function for cosmological parameters of H = 70, 
tt m = 0.3 and fi A = 0.7. 

L U: \(M) can be constrained from our data. Table 4 shows the sampled ?7-band luminos- 
ity within the common magnitude limit of My = —18.5 for all galaxies, for EL (star forming) 
galaxies and for NEL (quiescent) galaxies. All clusters have been sampled to at least this 
limit. We discuss possible contributions from fainter galaxies in App. C, but find that, 
unless the GLF shows a strong upturn at the faint end, we have sampled most of the ?7-band 
luminosity of the clusters. The [/-band luminosity is not corrected for radial sampling limits 
either, but represents the total emission from a cluster's galaxies within our survey area, 
which is approximately one harmonic radius. The table also gives the cluster mass, M 2 oo, 
taken from Reiprich & Bohringer (2002) and based on ROSAT and ASCA measurements of 
the intracluster gas profile. We estimate errors in the sampled £/-band luminosity with a 
Monte Carlo algorithm. 

Beijersbergen et al. (2002) have determined a [/-band GLF for the Coma Cluster from 
statistical background subtraction. Integrating their GLF to M v = — 19 gives a cumulative 
luminosity of L v « 369 x 10 39 erg s _1 A -1 . At a mass of M 20 o ~ 13.6 x 1O 14 M , this luminosity 
is lower than suggested by the extrapolation from the three clusters in our sample, but on the 
same order of magnitude. Given various systematic differences between these two surveys 

- most notably, the different survey areas and the use of statistical background subtraction 

- we refrain from imposing this data point as an additional constraint on our L l /(M 2 oo) 
relation. 

To derive an estimate for Lu t \ over a continuous range of cluster masses, we plot the total 
luminosities from Table 4 in Fig. 8 and fit an analytical expression to the three data points 
for the uncorrected cumulative luminosity. For similar analyses in the i^-band, Kochanek 
et al. (2003) get L K oc jVf L10±a09 , consistent with a linear relation. Girardi et al. (2000) 



Table 4. 



cluster 


mass [10 14 M ] 


Ljj sampled 


Ljj extrapolated 


Ljj extrapolated 
from EL/NEL 


LjJ L sampled 


Ly Bi sampled 


composite 
A496 
A754 
A85 


3.061°;" 

18.44^ ;o i 

7 61 +°' 79 
' - D1 -0.73 


533li 
3211H 

76lli 
5561*1 


761 + 80 
' Di -70 

458+tI 
79511" 


10481S 4 

1971+1408 

H56l^f 


8911° 
7611 

71±1| 
109111 


U7+f 9 
24llf 6 
682l 7 l 
454+H 



Notes: Luminosities in 10 39 erg s 1 . Errors are based on lcr uncertainties in analytical fits. Sampled luminosities are for 



galaxies with Mjj < —18.5, extrapolated Mjj < —10. Masses are from Reiprich & Bohringer (2002). 
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Fig. 8. — Integrated [/-band luminosity as a function of cluster mass. Filled circles show 
contribution from the sampled magnitude range (My < —18.5), empty circles show extrap- 
olation to Mjj < —10 using the Schechter fit to the composite GLF. Solid and dashed lines 
are various analytical fits, which are discussed in §4.4. 
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find a weak, but significant departure from linearity for the .Bj-band. Their estimate is 
oc M 8 — M ' 9 . Therefore, previous work suggests both linear and curved relations. To 
take this uncertainty into account, we attempt to fit both to the data in Fig. 8, and provide 
below the range of results from both approaches. 

A linear relation does not fit the three clusters in our sample well. The best fit slope in 
units of 10 24 erg s" 1 M" 1 is 57 ± 0.7, with X 2 = 95. 

For the curved relation, in order to avoid a singularity of the mass-to-light ratio at the 
origin, we do not fit a power law, but opt instead for the functional form 

L U:X = a 1 lg(a 2 M + 1). (4) 

This functional form is not theoretically motivated, but it reproduces the data points and 
provides a plausible interpolation between them, passes through the origin, and is differen- 
tiable there. The fit is excellent. With a 2 = 0.63 x 10~ 14 Mq 1 , the best-fit value, we find 
ai = 709ljj| * 10 39 erg s" 1 A" 1 . 

With these two fits, we obtain luminosity densities from clusters in the mass range from 
A496 (M 200 « 3 x 10 14 M ) to A754 (M 200 ~ 1.8 x 10 15 M ) from 1.7 to 2.8 x 10 36 erg s" 1 
A -1 Mpc~ 3 . The quoted range comprises the difference between the two analytical forms 
as well as the la errors associated with each. This luminosity density contributes to the 
[/-band background intensity by ^ = 4.3 — 6.4 x 10~ 12 erg s" 1 cm~ 2 sr _1 A -1 Gyr -1 , which 
corresponds to an increase of 0.1 — 0.3% relative to the current background per Gyr. The 
higher estimate is associated with the curved functional form. 

Estimating the cumulative contribution of clusters to the [/-band background would 
require models for the evolution of the cluster mass function, the star formation rate, and the 
spectral energy distribution of star forming galaxies at high redshifts. A detailed discussion 
exceeds the scope of this paper, and for this reason, we only give the rate of change in 
the ?7-band background due to present-day clusters here. However, even if star formation 
rates in clusters and their progenitors were higher by an order of magnitude in the past 
(approximately the difference between star formation in present-day field galaxies and the 
peak of star formation history around z«l (Madau, Pozzetti & Dickinson 1998; Steidel et al. 
1999)), the total cumulative contribution of clusters at z < 3 to today's [/-band background 
would still be only a few percent. 

Clusters do not contribute much to the mass density of the universe either. Integrating 
the mass function from Jenkins et al. (2001) shows that clusters in the mass range covered 
by our sample contribute ~ 2% of the critical density, or ~ 8% of the matter density of 
the universe. On the other hand, a significant fraction of the [/-band background may be 
contributed by galaxies in virialized systems of lower mass, by isolated field galaxies, or by 
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non-stellar sources. Integrating the mass function above shows that groups and clusters with 
M 2 oo > 10 12 M Q account for a mass density of Vt m = 0.21, and additional contributions could 
come from isolated galaxies and groups with M 200 < 10 12 M . 

To estimate the contribution from groups with M 20 o > 10 12 M , we extrapolate the 
mass-L[/ relations found above and integrate the luminosity density between M 20 o = 10 12 
Mq, the lower bound of the mass of the Local Group, and a variable upper mass cutoff. Fig. 9 
shows the resulting integrated luminosity density as a function of the upper integration limit. 
The total luminosity density from groups and clusters with M 2 oo > 10 12 M & approaches 
15.0 - 53.7 x 10 36 erg s" 1 A" 1 Mpc~ 3 . Groups of ~ 10 12 to 10 13 M contribute most to 
this emission. The contribution to the ?7-band background from the entire mass range is 
^ = (42.7 — 138) x 10~ 12 erg s _1 cm~ 2 sr _1 A -1 Gyr -1 , corresponding to an increase of 
1 — 6% of the current [/-band background over 1 Gyr, depending on the functional form of 
the fit and the value of the current ?7-band background used. 



5. CONCLUSIONS 

We have calculated galaxy luminosity functions (GLFs) from total [/-band magnitudes 
in three nearby clusters. Our analysis is based on a spectroscopic sample, providing clus- 
ter membership confirmation for each galaxy and avoiding the need to resort to statistical 
background subtraction. 

We have introduced a new variant of the maximum likelihood method for calculating 
luminosity functions. Conventional maximum likelihood methods use binned or analytic dis- 
tributions in just one or two variables (e.g., absolute magnitude and surface brightness) as 
an ansatz and require the algorithm to be customized and rerun every time the form of this 
ansatz is changed, making it inconvenient to analyze galaxy properties in higher-dimensional 
parameter spaces (e.g., as a function of absolute magnitude in different filter bands, surface 
brightness, environment, morphology, star formation, etc., simultaneously). Our Discrete 
Maximum Likelihood method, on the other hand, does not assume a specific form or dimen- 
sionality for the ansatz a priori, but assigns a statistical weighting factor to each sampled 
galaxy. This weighting factor contains all completeness and volume corrections. Therefore, 
the full range of observables of each galaxy can be retained throughout the calculation and 
subsequently analyzed without having to customize and rerun the algorithm with a new 
ansatz for every new analysis. The DML method is therefore ideal for analyzing modern sur- 
veys that gather a large number of photmetric, spectroscopic and morphological parameters 
about each object. Like other maximum liklihood estimators, the DML has the advantage 
of being unbiased by density inhomogeneities in the sample. 
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Fig. 9. — [/-band luminosity density from groups and clusters as a function of the upper mass 
limit M 20 o- Lower mass limit for the integration is 10 12 M . Mass function is from Jenkins 
et al. (2001). The upper line corresponds to the extrapolation using a curved functional 
form, as discussed in the text. The lower line corresponds to a linear relation between M 2 oo 
and L\j. Uncertainties (shaded regions) are based on the uncertainties in the normalization 
of Lu — M 20 o relation. 
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The results of our GLF analysis are summarized below: 

• The [/-band GLF in clusters down to M v < —17.5 M v + 2) can be described by a 
Schechter function with M v = -19.81 ± 0.27, a v = -1.09 ± 0.18, and ^ = 142 ± 5 
galaxies cluster -1 mag -1 . 

• We compare the [/-band and i?-band GLFs and find that, although the difference is 
too subtle to be reflected in the Schechter parameter a, the ratio of faint (M > M*) 
to bright (M < M*) galaxies is significantly larger in the [/-band. This indicates that 
cluster galaxies are bluer at fainter magnitudes. 

• For quiescent galaxies, we find best fit Schechter parameters of M v = — 19.771q;3q, 
au = — 0.97io!i8> an d 4>*u — 133 ± 6. The star forming GLF can be fit with a power 
law with a slope of a = 2.16io;]J. 

• If the Schechter fit to the overall GLF can be extrapolated to faint magnitudes {Mu < 
— 10), we have sampled ~ 85% of the total [/-band light from the clusters within the 
limit of Mu < —17.5. Quiescent galaxies dominate the [/-band flux at M v < —14. If 
there is a faint end upturn in the overall GLF (as suggested by extrapolating the star 
forming GLF) that continues past M v = —14, we cannot rule out that dwarf emission 
line galaxies dominate the total [/-band flux. The uncertainty is primarily in the shape 
of the star forming GLF. 

• Rich clusters in the mass range of our sample (~ 3 x 10 14 to 1.8 x 10 15 M ) account for 
a [/-band luminosity density of (1.71 ±0.02) to (2.46 ±0.33) x 10 36 erg s -1 A -1 Mpc -3 . 
This corresponds to an increase of the current [/-band background of 0.1 — 0.3% per 
Gyr. Additional contributions from galaxies outside our spatial survey limits could 
increase this estimate by a factor of 2-4. Galaxies beyond our faint magnitude limit 
could add another ~ 20% if the overall GLF can be extrapolated to Mu — 10. 

We thank the anonymous referee for a thorough reading of the text and valuable sugges- 
tions. We also thank Dennis Zaritsky for his comments, and Anthony Gonzalez for providing 
a numerical estimate of the cluster mass function. DC and AIZ acknowledge support from 
NSF grant # AST-0206084. This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronautics and Space Administration. 
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A. THE DISCRETE MAXIMUM LIKELIHOOD METHOD 

Our aim is to recover the galaxy luminosity function, i.e., the parent distribution from 
which those galaxies with both redshifts and [/-band photometry in our sample have been 
drawn. Reconstructing the parent distribution from a set of sampled galaxies requires two 
corrections. First, a volume correction is needed to account for the fact that bright galaxies 
are over-represented because they can be observed to greater distances and therefore over 
a larger volume. This is not a problem in volume-limited surveys (e.g., of only one cluster, 
where all galaxies are at approximately the same distance), but it is a concern whenever a 
GLF is derived from a magnitude-limited sample spanning a range in redshift. For example, 
in our case, a galaxy of a given absolute magnitude may be observable in one cluster, but not 
in the two other, more distant ones. Second, a completeness correction is necessary because 
not all galaxies that are photometrically detected in our i?-band master catalog have been 
sampled, i.e., have both spectroscopic and [/-band photometric information. There are 
various reasons why photometrically detected galaxies are not part of the sample: some 
faint galaxies were not targeted for spectroscopic observations, a small number were targeted 
unsuccessfully, and others were not detected on the [/-band images. The completeness of 
the sample relative to the i?-band master catalog varies primarily as a function of apparent 
magnitude and surface brightness. 

GLFs for individual clusters are usually calculated by simply weighting each galaxy by 
the inverse of its sampling probability, which is either known (in a strictly magnitude-limited 
survey) or can be recovered empirically (by comparing the number of spectroscopically sam- 
pled galaxies of a given magnitude to the total number of photometric detections). Com- 
posites are usually calculated by scaling and averaging the GLFs from several individual 
clusters. However, this method does not make optimal use of the information in a survey. 
The maximum likelihood solution for the parent distribution function is not only determined 
by the fact that a galaxy has been observed in one cluster, but also by the question whether 
or not it would have been observable in any of the other clusters in the sample. 

A method that is typically used for field galaxies, but that could be adapted for a 
sample of clusters, is the V/V max method, which weights each individual galaxy by the 
inverse of the volume over which it could have been observed in the survey. This method 
has two drawbacks: First, it implicitly assumes that the galaxy space density distribution 
is homogeneous. This is primarily a problem for the field GLF, but if applied to clusters, it 
also requires introducing an additional scaling to account for the fact that different clusters 
have different numbers of member galaxies. The second problem is that the method is 
impractical if the sample is not constrained by a fixed limiting magnitude, but by a fractional 
completeness that varies with apparent magnitude, which is the case for our cluster data. 
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The definition of the "volume" over which a galaxy would have been observable in such a 
survey is not straightforward. 

The most popular algorithms for calculating luminosity functions from a sample span- 
ning a range in redshift are based on maximum likelihood (ML) methods. ML algorithms 
assume a parent distribution characterized by a limited number of parameters. Taking the 
selection criteria (sample completeness and magnitude limits) of the survey into account, 
they then calculate the probability that the observed sample has been drawn from this 
assumed parent, and iteratively adjust the parameters to maximize this probability. ML 
methods are preferred because they are not dependent on assumptions about the redshift 
distribution of galaxies in space and are therefore unbiased by density inhomogeneities due 
to large scale structure. This advantage applies particularly to determinations of the field 
GLF, but ML algorithms are applicable to every galaxy sample regardless of the redshift 
distribution, including clusters. 

Several variants of the ML method exist that differ in the way they represent the parent 
distribution. The parametric maximum likelihood method (PML) by Sandage, Tamman & 
Yahil (1979) characterizes the parent as a Schechter function (Schechter 1976) with three 
free parameters: the shape parameters M* and a and the normalization constant 0*. The 
disadvantages of this method are that assumptions regarding the shape of the parent are 
required, and more complex, multi-dimensional analyses are impossible because the Schechter 
function is a function of only one variable, absolute magnitude. Cross & Driver (2002) have 
expanded the PML to include absolute magnitude and surface brightness. 

Efstathiou, Ellis & Peterson (1988), in the stepwise maximum likelihood (SWML) ap- 
proach, parametrize the galaxy parent distribution with binned distributions, a method that 
could, in principle, be expanded to an arbitrary number of dimensions (for an example of an 
adaptation of the SWML method to two-dimensional distributions, see Christlein & Zablud- 
off (2003)). However, in a sample of moderate size such as ours, with only a few hundred 
galaxies, the advantages of binning galaxies are negated by shot noise if the number of bins 
approaches the number of galaxies in the sample (although the effects of shot noise can be 
reduced by projecting the distribution back onto a subspace of lesser dimensionality). 

Both the PML and SWML methods suffer from the fact that a parametrized form of 
the GLF has to be assumed a priori, and that variables that are not explicitly represented 
by this form are discarded. They are thus not ideal for multi-dimensional analyses, i.e., 
for recovering the galaxy distribution function (GDF), which, in analogy to the univariate 
GLF, describes the abundance of galaxies as a function of multiple variables such as their 
luminosity in different filter bands, surface brightness, star formation indices, environment, 
morphology, etc. 
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In our case, the fact that we are modeling the sample completeness as a function of tur 
and fiR, while plotting the luminosity function as a function of Mjj, requires us to keep track 
of a minimum of three variables (My, Mr, /j,r) and to treat the galaxy distribution function 
as a function of these (note that in this context fi is measured in the galaxy rest frame). The 
PML method is clearly ruled out for our purposes because there is no generally accepted 
functional form to describe the distribution of galaxies in (Mr,i71u, ^r) space. 

While the SWML could, in principle, be applied to this problem, it is not very effi- 
cient, particularly with regard to future work aimed at expanding our analysis of the galaxy 
distribution function to higher-dimensional parameter spaces (including such variables as 
star formation, morphology, and local environment). It either requires modifications to the 
algorithm to change the ansatz for the galaxy distribution function every time a different 
cut through the distribution function is to be made, or the distribution has to be binned a 
priori in as many dimensions as photometric, spectroscopic and morphological parameters 
are available for each galaxy. Both approaches are computationally cumbersome, because 
the ansatz has to be hardwired into the algorithm and modified every time the dimensional- 
ity of the distribution function is changed. Furthermore, there is no advantage to binning a 
galaxy distribution in so many dimensions that many bins hold only a very small number of 
galaxies. In such cases, higher-order moments of the galaxy distribution within a single bin 
can become important, so that the mean properties of the bin do not accurately represent 
the properties of the individual galaxies contained therein. 

We have therefore developed a new variant of the maximum likelihood algorithm that 
does not require any assumptions about the functional form of the galaxy distribution prior 
to the calculation, nor does it require binning the galaxies or sacrificing any dimensions of 
the available parameter space in the sample. 

In our algorithm, the GDF is represented by the sampled galaxies themselves, and not 
by bins or functional parameters that only model certain aspects of the sampled galaxy 
population. Our approach can be regarded as an extreme application of the SWML method 
for infinitesimally small bins in parameter spaces of arbitrary dimensionality. We formally 
describe the GDF as a sum of weighted delta functions. It is nonzero at all coordinates in 
parameter space where our sampled galaxies lie. For this reason, we refer to it as the Discrete 
Maximum Likelihood method (DML) throughout this paper. The ansatz for the GDF is: 



where C is a normalization constant and x n is a parameter vector for galaxy n of arbitrary 
dimensionality. In our case, the parameter space will be (Mr, M v , /jr). The weighting factors 
u n are free parameters to be determined by the DML algorithm. We use the nomenclature 
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(p to indicate that this ansatz can be generalized to arguments x of any dimensionality. This 
ansatz is similar to that of the C method (Lynden-Bell 1971; Choloniewski 1986), but the 
procedure for solving for the free parameters, iv, is different, with our method retaining the 
benefits of ML estimators. 

The decisive advantage of this method over the SWML method is that, instead of ab- 
sorbing the individual galaxies into an ansatz for the GDF that represents only a limited 
number of galaxy properties, it associates a weighting factor with each individual galaxy. 
The free parameters remain tied to the individual galaxies, rather than to fixed grid points 
in a pre-selected subspace of parameter space. Therefore, it provides a completeness and 
volume correction, but still retains the full range of photometric, spectroscopic and morpho- 
logical properties of the galaxy for a subsequent analysis. Computationally, the algorithm 
is independent of the dimensionality of the distribution function that is to be analyzed, and 
only needs to be adjusted to account for the correct dependencies of the sampling fraction. 
This combines the advantages of the V/V max method with those of the ML methods. The 
algorithm is also very simple to implement computationally. 

The ansatz in Eq. Al is not of course a realistic physical representation of the par- 
ent distribution function. After calculating the weighting factors, the distribution can be 
smoothed or binned, or treated with other multi-variate analysis techniques. In contrast to 
conventional techniques, these procedures do not have to be hardwired into the algorithm a 
priori, but can be applied after the completeness and volume corrections have been carried 
out, allowing us to perform a multitude of statistical analyses on a data set without having 
to find a new maximum likelihood solution for each analysis. 

The derivation of an algorithm to solve for the free parameters u is analogous to the 
SWML method. We start with the probability that a galaxy i with parameters x^ would 
have been observed in the survey: 



where Fj describes the field of galaxy i and includes redshift, distance, Galactic extinction, 
position in the sky, and other information. f(x"i \ Fj) is therefore the probability that a 
galaxy with properties Xi would have been included in the sample if it was in the position of 
galaxy j. All indices refer to specific members of the sample of both spectroscopically and 
[/-band photometrically sampled galaxies, i.e., u>h is the weighting factor associated with 
galaxy h. The 5 function now allows us to consider the function only at discrete points and 



Pi = 




(A2) 
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thus solve the integral: 

Pt = jt^-et, • (A3) 

We now form the composite probability to have obtained the observed sample from the 
assumed distribution and take its logarithm: 

InC = ^ ln(y] u g 8(xj, x g )) + ^ lnf(xj \ Fj) - ^ In ^ u g f(x g \ Fi). (A4) 



Taking the derivative by u h , we find 
dlnC _ \^ S(xi] x h ) \- f{x~h\ Xj) 1 \- /04 I Xj) 

d^h ~ ^Y, g Ug&(Xi; X g ) ^ J2 g Ugf(Xg I *i) ~ ^ , E fl I *i) ' 



Setting this expression to yields an expression for u h that is suitable for an iterative 
solution: 

? e/^/c^iVo ) • (A6) 



This algorithm converges very quickly. To obtain an absolute normalization factor, we 
integrate the GDF down to M v = — 17. We repeat this integration for each galaxy in the 
sample, applying its visibility conditions Fj to predict how many galaxies of a given M\j and 
Hn would have been sampled, had they been at the position of this galaxy. Thus, we derive 
a prediction for the total number of galaxies with Mjj < —17 that we would expect to have 
sampled, and compare it to the actual number to derive a normalization factor: 



C = N sampled (£ N;^ I F t ))-\ (A7) 

i 3 

where N c i^ is the number of galaxies sampled in the cluster to which galaxy i belongs, 
N sampled is the total number of galaxies in both the spectroscopic and £7-band sample, and 
f(xj | Fi) is the probability that galaxy j would have been sampled in the survey field and 
at the redshift and coordinates of galaxy i. 

Once the weighting factors have been calculated, it is possible to apply techniques of 
multi-variate analysis and plot different cuts through the GDF over any desired variable, 
or to restrict the analysis to subsets of galaxies, without having to modify the algorithm 
or recalculate the weighting factors. In principle, GLFs in the ?7-band and the F-band 
can be plotted from the same output, using the same weighting factors and normalization 
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constants. Furthermore, a sample can be split (e.g., by emission line properties), and the 
two distributions can be plotted and compared separately; the weighting factors do not have 
to be recalculated. Therefore, correct relative normalization between the two subsamples is 
automatically guaranteed. 

Care has to be applied in cases where, due to survey design, the sampling fraction 
model does not fully account for all dependences of the sample completeness. In our study, 
spectroscopic targets have been selected by approximate bj magnitudes. This introduces a 
residual dependence of the sampling fraction / on bj — R colors that is not accounted for 
in our approach, which only models / as a function of (mR, /Xr). If a subset of galaxies is 
selected by parameters that are correlated with bj — R colors, such as cluster membership 
or emission line properties, biases can result, because the mean empirical sampling fraction 
calculated for galaxies of a given (m R , ji R ) may be systematically different from the true 
sample completeness for the selected subset of galaxies. We resolve this problem by applying 
systematic corrections specific to the individual subsamples (particularly the EL and NEL 
subsamples) to the sampling fraction, which requires us to recalculate the weighting factors 
for these subsets. These corrections are discussed in App. B. 

For a similar reason, when calculating the GDF for individual clusters out of a sample 
of multiple clusters, it is not legitimate to first calculate the weighting factors for the entire 
sample and then just plot the galaxies associated with one particular cluster, because the 
weighting factors determined that way would contain volume corrections that are appropriate 
to a sample spanning a range in redshift, but not to the sample of an individual cluster. 
For a given absolute magnitude, the sampling fraction is different from cluster to cluster. 
Therefore, the weighting factors have to be recalculated for each individual cluster in order to 
analyze individual cluster GDFs. This is not a drawback of the DML method, but common 
to all ML methods. 

When the recovered value of the GDF is exactly zero in any region of parameter space 
(for example, for very faint Mr), it is necessary to understand whether this is because the 
sampling fraction is so low that none of a potentially large population of galaxies have been 
sampled, or because there are no galaxies with these properties. Again, this is a problem 
common to all LF algorithms, but a careful treatment is particularly important in the DML, 
because it represents the GDF solely with sampled galaxies and therefore, by default, does 
not probe the sampling fraction outside of the sampled regions of parameter space (whereas 
other methods like the SWML will return an undefined result if the sampling fraction is 
zero). 

The sampling probability for a galaxy of any given absolute magnitude, surface bright- 
ness, etc., can be probed in the DML by introducing mock galaxies. The critical value for the 
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sampling probability can be defined in the following way: Consider the the GDF in two di- 
mensions, Mu and Mr, and define a bin in Mjj. The number of galaxies in the well-sampled 
parts of the galaxy distribution in this bin, which is represented by the reconstructed GLF, 
is Nq. Further, assume that there is a region of parameter space within this Mjj bin with 
a low sampling probability. The number of galaxies in this region is N, and the sampling 
probability for such a galaxy is /. The poorly sampled region is a problem only if a) the 
expected number of sampled galaxies in this region is less than 1, i.e. Nf < 1, and b) the 
number of galaxies that exist in this region is a substantial fraction p of the total number 
of galaxies in this Mu bin, i.e. N > p(N + N ) > pN . This yields the condition for poor 
sampling, pN < f^ 1 (note that this expression is independent of the size of the hypothetical 
poorly sampled region). The reconstructed luminosity function is therefore not affected by 
a systematic relative error of more than p as long as the sampling fraction / is greater than 
the critical value of (piVo) -1 everywhere within this bin. By introducing mock galaxies with 
the given Mu, a value of /j,r characteristic for galaxies with this Mu, and arbitrary Mr, this 
requirement can be checked easily and the effective limits of the survey in Mr determined. 

Once we have found the effective survey limits (for example, in Mr), it is necessary to 
ascertain whether it is likely that substantial fractions of the galaxy distribution lie outside 
these limits. For example, for any Mu, there are extreme values of Mr for which no galax- 
ies could have been sampled spectroscopically. However, there is a fairly tight correlation 
between M v and M R , and for bright Mu, it is extremely unlikely that any galaxies reside 
beyond the faint effective limit in Mr. In §3.4, we use the procedure outlined above to define 
the effective absolute magnitude limits in Mr and Mu for our sample. 



B. DETERMINATION OF THE SAMPLING FRACTION 

This section describes the empirical determination of the sampling fraction, f(x* h | Fj), 
which is the probability that we have both [/-band photometry and spectroscopic information 
for a given object with certain physical properties x n (e.g., absolute magnitude and surface 
brightness) if it is in a particular field Fj (characterized by redshift, Galactic foreground 
extinction, position on the sky). We determine the completeness of the ?7-band and spectro- 
scopic catalogs relative to the i?-band photometric catalog. We assume the i?-band catalog 
to be complete for itlr < 18 and fiR < 23.61 (the approximate equivalent of a per-pixel 
detection threshold of 1.5(7 on a typical frame). This assumption is conservative; in fact, we 
detect objects more than five magnitudes fainter, but do not include them in the analysis. 
We cannot detect objects at lower surface brightness, but at m R < 18, the distribution of 
photometric detections in the (tur, hr) plane is well separated from this limit, indicating 
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that few objects are likely to be lost to the surface brightness limit. 

Following Bayes's theorem (Bayes 1764), we factorize the sampling fraction — the prob- 
ability that, for a given object, redshift information (cz) and [/-band photometry (U) are 
available on condition that the object meets our sample selection criteria. These selection 
criteria, denoted by sel, are the criteria that we apply to our galaxy sample to select those 
galaxies that are to be included in our GLFs. This is primarily cluster membership. For 
some analyses, we also select galaxies by their emission line properties, and comment on this 
at the end of this section. 

The expression for the conditional sampling probability is then 

p(cz A U | sel) = p(cz A U A sel)fp(sel). (Bl) 

We now transform Eqn. Bl into a form that is suitable for numerical evaluation and 
factorize it into separate terms for the spectroscopic and the [/-band photometric complete- 
ness. There are two straightforward factorizations that allow us to separate terms connected 
to the spectroscopic and the ?7-band sampling fractions: 

p(cz A U | sel) = p(U | cz A sel) x p(cz \ sel). (B2) 

or 

p(cz A U | sel) = p{U) x picz A sel \ U) x p(cz \ sel)/p(cz A sel). (B3) 

The separation of terms for the spectroscopic and [/-band photometric sampling frac- 
tions is motivated by the opportunity to apply specific, systematic corrections to either. It 
is only these corrections that introduce any difference between the two methods, but the 
differences are statistically insignificant (< 1% even in the faintest magnitude bins). We use 
the first factorization because of its greater simplicity. 

Both the spectroscopic and ?7-band sampling fractions cannot be described analytically 
in our case. We thus have to determine them empirically by comparing the number counts 
of all detections in the i?-band photometric master catalog with those with redshifts and 
[/-band photometry. Because the i?-band photometric catalog is complete, we calculate the 
sampling fraction as a function of rriR and /ir. These two parameters alone do not allow for 
an unambiguous description of the sampling fraction; color and/or environment also play a 
role. We discuss below how we address this problem. The bin widths that we use to calculate 
the sampling fraction are Am = 0.75 and A/i = 0.25. We have experimentally found that 
these choices of bin widths reproduce stable results for the resulting GLFs. 

p[cz | sel) cannot be determined exactly, because information about cluster membership 
is only available for spectroscopically sampled objects. The assumption p(cz \ sel) picz) is 
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only justified if the sampling fraction is not strongly correlated with the selection criteria for 
a given (mjj,^). Unfortunately, the experimental design of our survey — target selection 
in bj, compilation of a master catalog in R — violates that assumption, because there 
are systematic differences in B — R colors between cluster and field galaxies, and therefore 
systematic differences in the spectroscopic sampling probability for a given (ma, fin). This 
color selection bias could be remedied if the sampling fraction model explicitly contained 
color terms, but it is impractical to estimate the sampling fraction as a function of more 
than two parameters with a sample the size of ours. 

Instead, we use an algorithm that systematically compensates for color selection bias 
and estimates p(cz | sel) under the assumption that the targeting probability is uncorrelated 
with the selection criteria for a given m hj . The probability that a galaxy of a given m R 
obeying the selection criteria has been targeted is 

pitarget \ sel) \ m R = (Nl rget ^ nonsel N + + N+ rgetySel N° + 2N°N^ nonsel + 
2N° ,N + ) + i(N° ,N + + N + N° + 27V°7V + + 

target, sel I 1 \\ target,nonsel ' target, sel ' target,nonsel 1 

27V° ,N + ) 2 — 8N°N + (2N° ,N + + TV ,N + + 

target, sel I \ target, sel target,nonsel 1 target,nonsel target,nonsel ' 

NLget,selNLget,sel)) 1/2 /^N°N + . (B4) 

Here, N denotes the number of galaxies that have been photometrically detected and thus 
included in the _R-band master catalog. N target is the number of galaxies targeted for spec- 
troscopic observations. The subscript sel indicates galaxies that are to be included in the 
GLF (e.g., that are cluster members). The superscript indicates that these quantities are 
evaluated at m^, and the superscript + means that they are evaluated at m R + A(B — R). 
A(B — R) is the difference in B-R colors between galaxies that have been selected (e.g., that 
are cluster members), and those that have not. We roughly estimate the m#-dependent dif- 
ferences in AB — R between cluster and field galaxies from the sampled galaxies themselves. 

In all cases, bj is not an accurate photometric magnitude, but rather the approximate 
magnitude that served as the basis for the target selection for the spectroscopic survey. 
Furthermore, for this derivation we assume N targettSe i ps N target N cz>se i/N cz , i.e., that the 
spectroscopic success rate for a given (m#, /x#) is not strongly correlated with the subsample 
selection criteria. Any error in this assumption is likely to be small, because the success 
rate is very high for the magnitude range analyzed. We then modify the sampling fraction 
fs{ m R, Hr) by the factor pitarget \ sel) /pitarget) \ mR . We generally apply these corrections 
to our results, except where we explicitly specify otherwise in order to estimate the size of 
the effect. 



Ultimately, we cut our galaxy sample not just by cluster membership, but also by 
emission line properties. The color selection bias described above applies to this case as well, 
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because emission line (EL) galaxies have systematically different (B — R) colors than non- 
emission line (NEL) galaxies. Therefore, the probability that a galaxy has been targeted 
for spectroscopic observations is systematically higher if it is an EL galaxy. We remedy 
this problem in the same way as described above, by assuming a A(B - i?) color difference 
between EL and NEL galaxies and applying our systematic correction to recover an estimate 
for p(cz | sel). 

In principle, a similar correction would be required to account for the fact that the 
completeness of the [/-band catalog has a residual dependence on U — R. Although Fig. 2 
shows our compensation for the incompleteness of the [/-band sample to work very well in 
reproducing the correct i?-band GLF, systematic errors arising from this problem are more 
likely to affect the [/-band GLF. To quantify the impact that these errors may have on our 
results, we pursue the following approach: 

We model the completeness of the ?7-band photometric sample, as 

p(U | cz A sel) = fu = fu(m R , ray, fx R ) 

Ft \ , ®fu i - \ 

~ fu{m R , fX R ) + g^( m u - m U\ R ) 

= fu(m R , hr) + §^^( m u ~ rhu\ R ). (B5) 

We determine from our default sampling fraction model and obtain ^ (which 
is of order unity) directly from the sample, (mu — ray \ R ) is the difference between the 
mjj for which we determine the sample completeness and the mean ray appropriate for the 
given niR. Therefore, this approach uses a first-order approximation to recover the sample 
completeness for a given ray and /j,r instead of m>R and /j, R (thus removing most of the bias, 
as the dependence on apparent magnitude is stronger than that on surface brightness). 

The resulting [/-band GLF is consistent with the [/-band GLF without these corrections 
(p x 2 = 0.994). Because the systematic effects are negligible and the above argument is only 
an approximation to assess the order of magnitude of the effect, we choose not to apply these 
corrections in our calculations. 



C. HOW MUCH LIGHT FROM FAINTER GALAXIES? 

Our photometric and spectroscopic samples are limited in apparent magnitude and sky 
coverage. Because of these constraints, we have sampled only a part of the total galaxy 
population in each of our three clusters. In this and the following section, we will estimate 
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what additional contributions in terms of numbers of galaxies and cumulative luminosity 
could come from cluster galaxies outside our photometric and spatial survey limits. 

Our composite GLF reaches to M v = —17.5. For a GLF with a rising faint end slope, 
such as the one in our sample, there are many more faint than bright galaxies. On the other 
hand, most of these are dwarfs that contribute only little to the total [/-band light from a 
cluster. How much of the total [/-band light have we sampled within our survey limit? To 
quantify this fraction, we fit a functional form to the sampled data and extrapolate it beyond 
the magnitude limit. As functional forms, we use our Schechter functions for the overall and 
NEL GLFs, and power laws for the EL GLFs. There are two possible approaches: We can fit 
a Schechter function to the overall GLF and extrapolate it to a very faint cutoff magnitude, 
or we can fit functional forms to the EL and NEL GLFs, extrapolate them individually, and 
add the results. 

These two approaches are generally not equivalent, because the sum of the power law EL 
and Schechter NEL GLFs is not necessarily a Schechter function. While a single Schechter 
function, as we have fitted to our overall GLF, assumes that the faint end is approximated 
by a single power law, the superposition of an EL and NEL GLF — assuming that there is 
no break or turnover in either GLF — usually shows a strong upturn at the faint end, just 
beyond our sampling limit, due to the increasing dominance of the power law EL GLF over 
the Schechter NEL GLF with its flat faint end slope. As we do not know whether the overall, 
EL, or NEL GLFs can be extrapolated all the way to the lower end of the galaxy magnitude 
range, or whether there is a break or turnover at a magnitude fainter than our survey limit, 
it is not clear which of these approaches best represents the real GLF of cluster galaxies. 

Claims of an upturn at the faint end of the cluster GLF have been made in the past on the 
basis of photometric surveys that employed statistical subtraction rather than spectroscopic 
membership confirmation for background decontamination (Valotto, Moore & Lambas (2001) 
and references therein). Recent spectroscopic surveys of the GLF in clusters (Christlein & 
Zabludoff 2003; de Propris et al. 2003) have not turned up evidence for such an upturn 
in the R— and 6j-bands, suggesting that previous claims may have been affected by biases 
inherent in statistical background subtraction (Valotto, Moore & Lambas 2001). Some recent 
photometric surveys of nearby clusters (Trentham & Hodgkin 2002) also indicate that the 
faint end slope is substantially flatter than our estimate for the EL GLF. Nevertheless, the 
question is not yet settled, because the spectroscopic surveys are generally shallower than 
purely photometric ones and the results from photometric surveys are still ambiguous. The 
uncertainty in the slope of the EL GLF in our present study is too large to predict at which 
magnitude an upturn might become significant. 

Because of this uncertainty, we pursue both extrapolation approaches and compare the 
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results below. To estimate our uncertainties, we perform the extrapolations for a range of 
Schechter parameters within the la error contours for the individual fits to the overall, EL, 
and NEL GLFs. 

Fig. 10 shows the results for the cumulative [/-band luminosity per cluster (averaged 
over the three clusters in our sample) as a function of the limit M v to which we integrate 
the GLF. We use the zero-point flux for the Johnson [/-band given by Colina, Bohlin & 
Castelli (1996) to convert absolute magnitudes to luminosities. The luminosity of a source 
with Mjj = —19 is given by L v 19 = 2.065 x 10 39 erg s _1 A -1 . The faintest known dwarf 
galaxies with recent star formation in the Local Group have My ps —10 (Pritchet & van 
den Bergh 1999). We adopt this value in the [/-band as the faintest magnitude to which we 
extrapolate the GLF. 

Judging from the overall GLF, galaxies fainter than our survey limit of Mjj = —17.5 
contribute very little to the total cluster [/-band luminosity. Therefore, we have sampled 
most of the cluster [/-band light in our survey. Extrapolating the GLF to Mu = —10, 
the fraction of light in the bright end (M v < —17.5) is 0.85+o;q5 (errors based on the la 
errors in a v and M v ). It is notable that most (83%) of the [/-band light at M v < —17.5 
is contributed by giant NEL galaxies, which make up ~ 80% in number. NEL galaxies 
dominate the [/-band emission from clusters at least down to Mu = —14. 

However, if we extrapolate the EL and NEL GLFs individually, the sum of the extrap- 
olated GLFs, for Mjj > —17.5, is not as flat as the fit to the overall GLF. The sum shows 
a strong upturn at fainter magnitudes due to the dominance of the power law EL GLF at 
faint magnitudes. In this case, populations of EL galaxies much fainter than Mu = —17.5 
contribute substantially to the total ?7-band luminosity of a cluster if there is no break in the 
EL GLF. In fact, for a faint end slope a < —2, the luminosity integral diverges. If we allow 
for the existence of a steep faint end upturn that arises from the extrapolation of the EL GLF 
to fainter magnitudes, our best estimate for the fraction of the light sampled at Mu < —17.5 
is 0.62+^38- In the worst case, we would have sampled only 23% of the total [/-band light 
from the clusters, and EL galaxies much fainter than our magnitude limit would constitute 
the dominant source of the EL flux. If the cluster GLF is indeed flatter than our EL GLF 
power law slope, as suggested by the results of Trentham & Hodgkin (2002), then the real 
fraction of the light that we have sampled is likely to be near the top end of this range. 

Are these results valid for any individual cluster? The numbers above apply to an 
"average" of the three clusters in our sample. However, we have demonstrated in §3.1 that, 
even if the GLFs in our clusters are consistent with each other in shape, they differ in 
normalization. To estimate the cumulative Lu for any individual cluster, we therefore have 
to consider them separately. We use the composite GLF to constrain the shape parameters 
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Fig. 10. — Integrated [/-band luminosity for emission line (EL), non-emission line (NEL) and 
overall GLFs per cluster as a function of faint integration limit. The solid bold line shows 
the cumulative [/-band luminosity of all galaxies, using the Schechter fit for the overall GLF 
for extrapolation. The long-dashed line shows the same for NEL galaxies. The dotted line 
shows EL galaxies. The upper, short-dashed line gives the sum of the EL and NEL integrals. 
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for the extrapolation under the assumption (which is consistent with the data) that the same 
fit describes all three clusters. We then renormalize the Schechter function to reproduce the 
total luminosities at My < —18.5 in each individual cluster. Table 4 shows the luminosities 
for the composite GLF and the individual clusters. It gives both the total luminosities at 
the bright end {My < —18.5) and the extrapolated luminosities to a limit of My = —10. 
The table shows that the systematic uncertainties about the shape of the faint end GLF, 
particularly whether there is an upturn or not, are larger than the random uncertainties 
associated with either extrapolation method. 

Note that we have adopted a limit of My = —18.5 instead of the My = —17.5 that we 
regard as our confidence limit for the composite GLF. This is because, for the same reasons 
laid out in §3.4, the more distant clusters, A754 and A85, are only sampled to M R & —19.0. 
For My > —18.5, significant numbers of galaxies could be lost at fainter M R . However, we 
still use the full extent of the composite GLF (to My = —17.5) to constrain the Schechter 
parameters for the extrapolation. 

D. HOW MUCH LIGHT FROM LARGER RADII? 

What is the effect of the limited spatial coverage of our survey? We calculate galaxy 
number and luminosity profiles as a function of radius (normalized so that the harmonic 
radius is exactly 1 in each cluster). We then fit these profiles with isothermal beta models 
(Cavaliere & Fusco-Femiano 1976, 1978) projected onto the radial dimension, and extrapo- 
late the models to estimate the total number of galaxies and total U -band luminosity out to 
4.1 harmonic radii (~ 5 Mpc in the case of A496, roughly the infall radius). This extrapo- 
lation is justified by recent results from Gomez et al. (2003), who show that star formation 
rates typically exhibit a break around 3-4 virial radii, comparable to the outer limit of our 
extrapolation and to the edge of the cluster. 

We estimate that our survey samples 0.50^Q og of the total [/-band luminosity of a 
cluster from galaxies with My < —17.5 within the survey region of A496. The fraction of 
galaxies sampled in this region is 0.80^g'g. This result indicates that the ?7-band luminosity 
profile is more extended than the galaxy number count profile, presumably due to radial 
color gradients in the cluster. Within 0.9 harmonic radii, our survey samples more than 
three quarters of the galaxies in our clusters, and about 1/2 of the f/-band luminosity from 
galaxies with My = —17.5. 
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Fig. 11. — Integrated £/-band luminosity as a function of the cluster sampling radius in 
units of the harmonic radius for each cluster. The dotted line shows the measured profile 
out to a limit of r = 0.97\,j r (denoted by a thin vertical line). The dark solid line shows the 
extrapolation from the best fit isothermal-/? model. Errors are based on the la uncertainties 
in the best fit. The extrapolated profile is normalized to yield the correct sampled luminosity 
within r = 0.9r vir . Within this limit, we sample ~50% of the [/-band light from galaxies 
with M v < -17. 
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